Resolution of Single Spin-Flips of a Single Proton 
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The spin magnetic moment of a single proton in a cryogenic Penning trap was coupled to the 
particle's axial motion with a superimposed magnetic bottle. Jumps in the oscillation frequency 
indicate spin-flips and were identified using a Bayesian analysis. 
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Recent dramatic advances in quantum control of a sin- 
gle isolated nucleus have opened the way for direct pre- 
cision measurements of the proton and antiproton mag- 
netic moments. The present most precise value for the 
proton magnetic moment comes from measurements of 
the hyperfine splitting in atomic hydrogen J^. Bound- 
state corrections have to be included to extract the mag- 
netic moment of the free proton with a precision of 8.2 
parts in 10^ The antiproton magnetic moment has 
been determined with a precision of 2 parts in lO'^ [3 
from super-hyperfine spectroscopy of antiprotonic helium 
[4] and from fine-structure spectroscopy of antiprotonic 
lead [g. 

A direct precision measurement with just one isolated 
particle in a Penning trap has the potential to improve 
the precision of the value of the proton magnetic moment 
/ip by more than one order of magnitude. In addition, 
there would be no need for theoretical corrections. Thus 
a direct measurement could be used for consistency tests 
with previous measurements based on hyperfine split- 
tings. In the case of the antiproton the potential for 
improvement is more than six orders of magnitude. This 
would enable another stringent test of the symmetry be- 
tween matter and antimatter in the baryon sector [5]. 
The principle of a direct measurement of fj,p is to deter- 
mine the Larmor (spin-precession) frequency vl and the 
cyclotron frequency Vc of a proton in a magnetic field. 
The frequency ratio vl/^c = fJ-p/li-N gives the proton 
magnetic moment fj,p in terms of the nuclear magneton 
UN- The cyclotron frequency i>c of a proton in a Pen- 
ning trap can readily be measured by preparing a double- 
dressed state [7] and applying the Brown-Gabrielse in- 
variance theorem [8]. The Larmor frequency vl can be 
determined by driving spin-flips and measuring the tran- 
sition probability as a function of the drive frequency. 
To this end, an inhomogeneous magnetic field, a "mag- 
netic bottle," is used which couples the spin magnetic 
moment to the axial motion of the proton. Using this 
"continuous Stern-Gerlach effect" [21 HO], jumps in the 



axial oscillation frequency indicate spin- flips. The chal- 
lenge is to detect these spin-flips on a background of axial 
frequency changes which result from tiny changes of the 
motional angular momentum of the proton in the trap. 
Recently we reported on the statistical detection of 
spin-flips of a single proton in an inhomogeneous mag- 
netic field [TT]. We have used this method to obtain 
8.9 • 10~^ [12], while a similar experiment by another 
group achieved 2.5 • 10~^ [13]. Both measurements are 
limited by the inhomogeneity of the magnetic field. The 
precision can be boosted by several orders of magnitude 
using the double- Penning trap technique [14 . In this el- 
egant method the measurement of i^c and the excitation 
of spin-flips at happens in a first Penning trap with 
a homogeneous magnetic field. Spin-state detection is 
carried out in a second Penning trap with a magnetic 
bottle. This method has been used to measure the mag- 
netic moment of the electron bound in ^*Si^'^+ with a 
relative precision of 5-10~^° [T3j. The double-Penning 
trap technique requires that single spin-flips can be re- 
solved, which so far was not possible with nuclear spins. 
In this Letter, we present the first detection of single 
spin-flips of a single proton. Noise-driven random transi- 
tions between the cyclotron quantum states cause a back- 
ground of frequency fluctuations. The characterization 
of these frequency fluctuations enables a novel spin state 
analysis method for Penning trap experiments based on 
a Bayesian formalism. 

Our apparatus consists of a cryogenic double-Penning 
trap, mounted in a superconducting magnet and cooled 
by a liquid helium cryostat. Both Penning traps are 
shown in Fig. [T] and have five electrodes in compensated 
and orthogonal design [12l US]. The precision trap is 
placed in the central homogeneous volume of the 1.89 T 
magnetic field. The analysis trap has a ring electrode 
made of ferromagnetic Co/Fe material, which shapes the 
magnetic field to a so-called magnetic bottle 
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FIG. 1. Schematic of the experiment. Two cylindrical Pen- 
ning traps are connected by transport electrodes. Supercon- 
ducting toroidal coils and cryogenic low-noise preamplifiers 
are used for the detection of the axial motion. The signals 
are analyzed by a fast Fourier transform, FFT. The central 
ring electrode of the analysis trap (black) is made of a ferro- 
magnetic Co/Fe alloy. The magnetic field along the z axis is 
indicated in the lower graph. Voltages are applied to the elec- 
trodes using low-pass filters (not shown). For further details 
see text. 



with B2 = 2.97(10) • lO^T/m^. The proton can be 
moved between both traps using transport electrodes. 
The whole electrode stack is placed in a sealed vacuum 
chamber and cryopumping is utilized. Collisions with 
residual gas are negligible and single particles can be 
trapped for months. 

In the Penning trap the proton has three eigenmotions, 
the axial motion with frequency v^., and two radial mo- 
tions: the magnetron and the modified cyclotron motion, 
at frequency v- and respectively. The eigenfrequen- 
cies are = 623 kHz, = 8 kHz, j/+ ^ 28.9 MHz in 
the precision trap, and Vz = 742 kHz, 1/^ = 15 kHz, 1/^ = 
17.9 MHz in the analysis trap. The motion of the pro- 
ton induces image currents in the trap electrodes. These 
currents in the fA range are measured by connecting a 
superconducting inductance L to the trap, which forms 
a resonant circuit together with the parasitic trap capac- 
itance. The resonant circuit has a quality factor Q and 
acts as an effective parallel resistance Rp = 2t:i' ^.^sQL on 
its resonance frequency Vres- The signals are amplified 
and a fast Fourier transform (FFT) is performed to access 
the particle's eigenfrequency. Both traps are connected 
to detection circuits for the axial motion. The detection 
system for the precision trap has a quality factor of 12500, 
resulting in Rp = 130 Mfi and a signal-to-noise ratio of 
S/N — 20 dB. Here the noise is given by TV = m„, with 
Un the voltage noise density of the preamplifier, and the 
signal by 5*^ = 4:kBTRpl\va'^ + i'^R'^a^ ^-u^, with i„ the 
current noise density of the preamplifier, l^u the band- 
width and a the decoupling of the inductor and pream- 
plifier [17] . The resonator for the analysis trap has a 
quality factor of 9500 with Rp = 85 MTL at S/N = 25 dB. 



At Vres = the trapped proton shorts the thermal noise 
of the detector. This causes a dip in the noise spectrum 
of the detector from which v^, can be determined jI8) . 
Additionally, one cyclotron damping coil is connected to 
a split electrode of the precision trap . 
The proton spin magnetic moment [ip is coupled to the 
axial motion by the magnetic bottle in the analysis trap. 
This causes an additional potential term ^mag = — 
which adds to the electric quadrupolar potential $^ of 
the Penning trap. \l is the sum of the magnetic moments 
/!+ and /i_ due to the modified cyclotron- and magnetron 
motion, respectively, and the spin magnetic moment \lp. 
The total axial potential energy of the single particle is 
-Epot = Ip^E ~ M-B and the axial frequency is 
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where C2 is a geometry parameter of the trap and Vq the 
trapping potential. Changes in fipBo, and \E- \ shift 
the axial frequency by 
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In our analysis trap a proton spin-flip causes an axial 
frequency jump of Ai'z,sf = 171 mHz at 742 kHz. The 
strong magnetic bottle, however, also makes the axial 
frequency extremely sensitive to the energy in the radial 
modes. A change of only 4 parts in lO'' of the thermal 
cyclotron energy of 360 ^eV (4.2 K) causes the axial fre- 
quency to change by the same amount as a proton spin- 
flip. 

Therefore, for the detection of individual spin-flips the 
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FIG. 2. (a) Cyclotron transition rate as function of the cy- 
clotron energy, (b) Measurement of the axial frequency sta- 
bility as a function of averaging time. The dashed line which 
decreases is due to detector white noise averaging. The dot- 
ted line comes from a calculation with white noise driving the 
cyclotron mode (see text). The lower horizontal line indicates 
the best stability of Sopt = 55 mHz achieved. The upper hor- 
izontal line indicates the size of the frequency jump due to a 
spin-flip of 171 mHz, which is about 3Hopt. 

stability of the radial energy is crucial. Transitions be- 
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tween the cyclotron quantum states n+ cause axial fre- 
quency jumps of ^Vz,+ = i63 mHz. As a result the axial 
frequency fluctuates with "E^^y^ — {Sn+/St)TAi'^^^, where 
T is the FFT averaging time. The cyclotron transition 
rate can be calculated using Fermi's golden rule 



~5t 



27r , 



(4) 



Here p{E+) is the density of states of the 1-dimensional 
harmonic oscillator and 27rA+ = 6-82 (z^) /rrip is the 
line- width of the cyclotron resonance due to the cou- 
pling of the axial mode to the thermal bath of the detec- 
tion system, (z^) is the expectation value of the square 
of the axial amplitude. 
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is the electric dipole cyclotron transition matrix element 
with Eq the spectral noise density of a spurious elec- 
tric field. Note that Eq. Q gives a cyclotron transi- 
tion rate which is proportional to the cyclotron energy, 
Sn^/St (X n+. 

In the experiment E^ was calibrated by determining 
Vz + lS.Vz{Ej^) with the method described in |2l]. The 
proton's cyclotron mode was thermalized with the damp- 
ing coil in the precision trap several times, which results 
in Boltzmann distributed frequency shifts l^Vz{E^) in 
the analysis trap. E+ = is then identified by Ai^^ — 0. 
After the calibration Scyc was measured for different cy- 
clotron energies and the transition rate 5nj^ /St was de- 
termined as a function of E^ . The measurement data are 
shown in Fig. and are consistent with a linear rela- 
tion as predicted above. {5n+ / 5t) / E^ = 0.35meV~^s~^ 
was obtained, which corresponds to a noise drive of 
Eq =7nV-m"i-Hz"i/2 ^o obtain S^yc < Ai'z.sf/'i at 
a typical FFT averaging time of 200 s a proton with a 
cyclotron energy of less than 10 /leV has to be selected. 
As a measure for the axial frequency fluctuation 
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FIG. 3. Comparison of the threshold method and Bayes anal- 
ysis, (a) shows contour lines of the fidelity for a spin-flip 
probability of 50% as a function of the white noise (t„ and ran- 
dom walk ffu, contribution; solid hues: Bayes analysis, dashed 
lines: threshold method, (b) shows the fideUty as a function 
of the spin-flip probability at = (t„ = 39 mHz. 

we define S as the standard deviation of the dif- 



ference between two subsequent axial frequency mea- 
s urements a[T) = Vz{t) - Vz(t + T), S(r) = 
^(7V - l)-i Y.j^{a{T) - a{T)Y. E is invariant under 
drifts and, in the case of {ct{T)) = 0, has the same charac- 
teristics as the common Allan deviation [22 . InFig.[2];b) 
the variation of S with the measuring time T is shown. 
For short measuring times S decreases with 1 / -s/T due to 
averaging of the detector white noise as indicated by the 
dashed line. At measuring times above 300 s 5 increases 
with VT, which is due to a random walk caused by fluc- 
tuations of the cyclotron energy. The dotted line is the 
result of a calculation with a transition rate for cyclotron 
quantum jumps of 6n+/5t = 0.002 s~^. 
The lowest axial frequency fluctuation achieved is Eopt = 
55 mHz, which corresponds to a reduction of the fre- 
quency fluctuations due to the random walk and white 
noise by about 90 % compared to our previous work [T^] . 
This improvement is due to the detection system with a 
higher quality factor reducing the white noise contribu- 
tion. By an increase of the effective electrode distance 
a stronger decoupling [THj is possible. Because of the 
smaller line- width of the noise dip and the higher signal- 
to-noise ratio of the detection system the axial frequency 
is measured faster and more precisely. For a reduction 
of a cyclotron noise drive a split electrode of the analy- 
sis trap, initially used for coupling of the axial mode to 
radial modes, was replaced. 

Spin state analysis can be done with a simple method 
assigning each axial frequency jump above a given thresh- 
old to a spin-flip. We developed an alternate and more 
advanced method, which is based on probability theory 
and uses an update probability theorem, Bayes rule |23| . 
In a series of measurements, each axial frequency fi is 
modeled as the sum of an accumulated random walk 
Wi = ^Wi, white noise Ui and a frequency jump due 
to a spin-flip with fi = Wi + rii ± 1/2/S.Vz.sf ■ The distri- 
butions of Wi and rii can be described by normal distri- 
butions N{^] a) with mean /i and standard deviations 
a^i and cr„, respectively. The conditional probability 
Pa,w = P{(^i,Wi\fi,fi-i,--) is defined, which is the 
probability of being in a spin state = {f i, ii} and hav- 
ing an accumulated random walk Wi given all frequency 
information fi, fi-i, .... This definition allows the appli- 
cation of Bayes's rule relating the a-posteriori probability 
IV a-priori probability P (ofi, Wi \ /i-i, ...), 

P W„ ...) P (a„ W^,|/._i, ...) 
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with P{fi\fi^iT--) the normalization found by 
!w '^aPaW ~ 1- frequency correlation in- 

troduced by the random walk the conditional probability 
of measuring frequency fi after a measurement of fi-i 
cannot be regarded independently, P{fi\fi-{) ^ P{fi)- 
In our definition ai and Wi completely specify the 
distribution of fi and the first factor in the numerator 
can be simplified, P {f^\a^,Wi, f^-i, ...) = P{f^\a,,Wi). 
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FIG. 4. Observation of single spin-flips with a single proton. At the top a series of axial frequency measurements is shown. 
A resonant spin-flip drive was applied between the crossed data points, indicated by thick arrows. In between an off-resonant 
spin-flip drive was turned on indicated by the thin arrows. Otherwise the drive was turned off. Large axial frequency jumps 
after a resonant drive are due to single proton spin-flips. At the bottom the result of the Bayesian analysis is shown. 



Hence the first factor is the probability density for /i, 
given the state as well as random walk information at 
time i. It depends on the white noise contribution, 

P{Ma,,W,) = N{f, -W,± l/2Ai.,,,/;a„), (7) 

with + for ai —^i and — for respectively. The 

second factor in the numerator of Eq. [6] depends on the 
previous state through the previous frequency measure- 
ments /i-i,/i-2, ■•■ • Integrating over all possible states 
and random walks gives 

P ( (8) 




The first factor describes the probability density without 
the knowledge of the latest frequency measurement. It 
introduces an update of the state probabilities describing 
the evolution of the state probability from measurement 
j — 1 to i, 

P ( a„W^,|a,_i,W,'_i, /,_!,...) 
= F(a„PF,|a,„i,VF/_i) 

= {l-Psf)N{W,-WU-a^) fora, (9) 
= PssN{W^ - WU-a^) for a, ^ (10) 

with the spin transition probability psf. The second 
factor in Eq. |8| -P^~^^ simply is the state and random 
walk probability at time The final state probability 

at time i, P^, is obtained by usage of the marginalization 
rule = . Input parameters are the spin- flip 

probability Psf^ the random walk contribution and 
the white noise contribution (T„. 

To investigate the quality of the spin state analysis we 
define the fidelity as the fraction of correctly identified 
spin states in a series of frequency measurements. It was 
found numerically that the best fidelity in the threshold 
method is obtained for a threshold of Ai^^ ,j//2. 



Fig. [sja) shows the fidelity achieved with both, the 
threshold method, and the Bayes algorithm using 
simulated data. In the case of a pure random walk 
both methods give the same fidelity. If white noise 
contributes significantly the Bayes algorithm has better 
fidelity than the threshold method. In Fig. [sj^b) the 
fidelity is shown as a function of the spin-flip probability 
using again simulated data. The random walk and white 
noise contribution were set to values which correspond 
to the optimal frequency stability ((7^, = (T„) in the 
experiment, Sopt = y^cr^, + cr^ = SSniHz. The Bayes 
method clearly is superior to the threshold method and 
has a fidelity of 88 % at 50 % spin-flip probability. 
The experimental observation of single spin-flips is 
shown in Fig. |4] The upper trace shows a time-series 
of axial frequency measurements. Thin vertical arrows 
indicate times when an off-resonant radio-frequency was 
applied to the spin- flip coil (see Fig.[T]). No extraordinary 
frequency jumps are seen in the data at these times. 
This is an important background test which shows, that 
the spin-flip drive does not afl'ect the cyclotron motion. 
Thick vertical arrows indicate times when a resonant 
drive was applied. Several axial frequency jumps whose 
size corresponds to spin-flips can clearly be seen in the 
raw data exactly at these times. One example of a 
spin- flip down and later up is at 240 min and 260 min. 
The lower trace in Fig. |4] is the result of the Bayesian 
spin-state analysis which is initialized with uncertainty 
(50% spin down probability) at f = 0. Note, that the 
Bayesian analysis is causal |24j . i.e. at each point in time 
it uses only present and prior axial frequency data. The 
Bayesian analysis nicely confirms the spin-flips visible 
in the raw data and provides a consistent picture of the 
time-evolution of the proton spin-state projection. 
In conclusion we observed single spin-flips of a single 
proton for the first time. This enables the application of 
the double Penning-trap method to measure magnetic 
moments of both the proton and the antiproton with 
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10~® precision, or better. 
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